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ABSTRACT 

The X-ray binary pulsar GX 1-1-4 was observed with Ginga every year from 
1987 to 1991, and with ASGA in 1994. During the Ginga observations, GX 1-1-4 
was in the steady spindown phase, though the X-ray flux was not steady. 
Assuming a distance of 10 kpc, the absorption-corrected X-ray luminosity 
decreased down to L 2 _ 2 Q ~ ^ in 1991, after the peak 

activity of T2_20 keV ~ ^ in 1989. On the other hand, the 

absorption column density showed a drastic increase over the Ginga observation 
series. It was less than 10^^ cm“^ at the beginning of the series, and it reached 
a maximum of (1.4 ± 0.2) x 10^"^ cm“^ in 1991, indicating a rapid accumulation 
of matter in the vicinity of the source. The center energy and equivalent width 
of the iron line were consistent with emission by isotropically distributed cold 
matter. The ASGA observation was performed on 1994 September 15, a month 
before the transition into a spinup phase. The source brightened again to 
-^2-20 keV ~ The absorption column density was observed to 

decrease for the first time to (2.08 ± 0.02) x 10^^ cm“^. The ionization degree of 
iron in the absorbing matter was determined to be Fe l-Fe iv using the ratio of 
the line-center energy to the absorption-edge energy. The low ionization degree 
is consistent with an absorbing matter distribution extending ~ 10^^ cm from 
the source. We compared the results with optical observations, and found the 
optical data also supports the picture. Based on the geometrical model, possible 
causes of the bimodal behavior of the source are discussed. 

Subject headings: binaries: symbiotic — pulsars: individual (GX 1-1-4) — 
X-rays: stars 
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1. Introduction 

The X-ray pulsar GX 1-1-4 is the most enigmatic of all binary X-ray pulsars. It is 
identihed with the symbiotic system V2116 Oph, which contains a magnetized neutron star 
accreting from the M6III companion (Davidsen, Malina & Bowyer 1977). GX 1-1-4 is a slow 
pulsar with a spin period of about 2 minutes. It has a hard energy spectrum, and hence has 
been a good target for both soft and hard X-ray observations. 

In the 1970’s, GX 1-1-4 was a bright X-ray sonrce with a typical 2-6 keV flnx of 100 
mGrab (McGlintock & Leventhal 1989 and references therein). (Hereafter, the bandpass 
of flux measurement is 2-6 keV, if not specihed.) It showed the fastest spinnp among the 
binary X-ray pulsars with a mean rate of P/P = —2.7 x 10“^ yr“^ (Eisner et al. 1985). 

In the early 1980’s, the sonrce entered into a low state. EXOSAT observations in 1983 
failed to detect the sonrce and set an npper limit of 0.5 mGrab (McGlintock & Leventhal 
1989 and references therein). Ginga detected GX 1-1-4 in 1987 at a flnx level of 2 mGrab 
and the source was found to be spinning down (Makishima et al. 1988). Snbseqnent 
Ginga and HEXE observations showed that the spindown rate was remarkably constant at 
P/P = 1.3 X 10“^ yr“^ in spite of a gradual increase by an order of magnitnde in the X-ray 
flux (Sakao et al. 1990; Mony et al. 1991). The clear contrast of GX 1-1-4 in the 1970’s and 
in the 1980’s indicates the close connection between the high/low transition in Inminosity 
and the spin reversal. 

The spindown rate was observed to decrease a little in 1990 and 1991, althongh 
snbseqnently a linear trend of spindown was recovered with a slightly increased rate 
(Mandron et al. 1994). A large flare in the hard X-ray band was detected by GGRO in 1993 
September (Stanbert et al. 1995).and the sonrce brightened np to ~80 mGrab (40-100 keV) 
for a few weeks. A gradnal increase of the spindown rate had been observed for abont 6 
months before the flare (Finger et al. 1993), bnt the rate was constant dnring the flare at 
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P/P = 3.2 X 10 ^ yr ^ (Staubert et al. 1995). 

The spin period history of GX 1+4 has been almost continuously monitored by 
BATSE since 1991 (Chakrabarty et al. 1997). It revealed a short episode of spinup in 1994 
(Chakrabarty et al. 1994), which was accompanied by a large increase in the hard X-ray 
flux. 

The 20-100 keV phase-averaged pulsed flux reached ~100 mCrab in 1994 November. 
The spinup phase was terminated in 1995 March (Chakrabarty et al. 1995) and the 
spindown resumed at a rate of P/P = 4.8 x 10“^ yr”^- 

The origin of this enigmatic behavior in GX 1+4 is not yet understood. Even the basic 
parameters of the system, such as the orbital period and the surface magnetic held of the 
neutron star, are not known. A binary period of 304 day was suggested by Cutler et al. 
(1986) based on the huctuation of the spinup rate. However, optical monitoring of the Hq 
emission line indicates that the period may be signihcantly greater than 304 days (Sood et 
al. 1995). If the spin behavior is analyzed in the framework of the accretion torque model 
by Ghosh and Lamb (1979a, 1979b) and Wang (1987), a very large magnetic held on the 
neutron star surface, rsj 10^^ G, is required (Makishima et al. 1988). Dotani et al. (1989) 
showed that even a moderate magnetic held of a few xlO^^ G can explain the spindown rate, 
if a retrograde accretion disk were formed around the neutron star. An energy-dependent 
pulse prohle in hard X-rays is used to construct a model based on two-photon emission 
as the dominant continuum source (Greenhill et al. 1993). The retrograde disk model is 
supported by the BATSE observations (Nelson et al. 1997), although some other possibilities 
such as the transition between Keplerian how and advective how in accretion geometry are 
also suggested (Yi, Wheeler & Vishniac 1997). 

In this paper we report results obtained from the unpublished Ginga and new ASCA 


data. 
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2. Observations 

The data we analyzed in this paper were obtained with Ginga (Makino & the Astro-C 
Team 1987) and ASCA (Tanaka et al. 1994). The main instrnments aboard Ginga are 
the Large Area Gonnters (LAG; Turner et al. 1989), which are non-imaging, conventional 
proportional counters. The LAG has an effective area of 4000 cm^ and covers an energy 
range of 1-37 keV or 2-60 keV depending on the high voltage setting It has the relatively 
poor energy resolution of ~18 % at 6 keV. The field of view of the LAG is restricted by 
collimators to 1° x 2° (FWHM). 

Ginga observations of GX 1-1-4 were made on 5 occasions. The first 3 sets of data were 
already analyzed and the results were previously published (Makishima et al. 1988; Dotani 
et al. 1989; Sakao et al. 1990). We used the 4th and 5th sets of data for the present analysis. 
The fourth observations were carried out on 1990 September 4-5, as a part of galactic ridge 
observations, and the 5th ones were on 1991 September 16. In the 5th observation, the high 
voltage setting of the LAG was changed to cover 2-60 keV, while in the 4th observation the 
nominal setting of 1-37 keV was adopted. Both sets of data were acquired in MPGl mode 
with 48 energy channels. In this data mode the temporal resolution was 0.5, 4, and 16 s for 
high, medium, and low bit rate, respectively. Unfortunately, most of the 5th observation 
was recorded in low bit rate mode, and this limited the temporal resolution in later analysis. 

ASGA carries four identical grazing-incidence X-ray telescopes (XRT), and each XRT 
has an imaging spectrometer at its focal plane. The XRT utilize nested thin-foil conical 
reflectors to achieve large effective area (1200 cm^ at 1 keV; Serlemitsos et al. 1995). A point 
source image has a typical size of 3' in half power diameter. The focal plane detectors are 
two solid-state spectrometers (SIS; Burke et al. 1994) and two gas imaging spectrometers 
(GIS; Makishima et al. 1995). The SIS consists of 4 GGD chips optimized for X-ray photon 
detection. It covers an energy range of 0.5-10 keV and a single chip covers a IT x 11' sky 
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field. The SIS achieves high energy resolution of 2 % (FWHM) at 6 keV. The GIS is a 
gas scintillation proportional counter with imaging capability. It covers an energy range of 
0.7-10 keV and the held of view of 25' in radius. The GIS has an energy resolution of 7.8 % 
(FWHM) and a spatial resolution of 0.5 mm (FWHM) both at 6 keV. 

ASGA observations of GX 1+4 were carried out in 1994 September 14-15. The 
observation mode of SIS was 1-GGD faint mode, and that of GIS was PH mode with 
standard bit assignment. Thus the time resolution of SIS was 4 s and that of GIS was 62.5 
ms and 0.5 s for high and medium bit rate, respectively. The journal of the observations are 
summarized in Table 1. 


EDITOR: PLAGE TABLE | HERE. 

3. Analysis and Results 
3.1. Ginga data 

Raw data were hrst screened according to the good time intervals, which exclude South 
Atlantic Anomaly passages. Earth occultation of the source, and the periods of low cut-off 
rigidity (<6 GeV/c). Light curves of the data were then checked to remove parts of the 
data contaminated by particle events manually. These screened data were used for the 
subsequent analysis. 


3.1.1. Timing analysis 

The pulse period of GX 1+4 was determined using an epoch folding technique. We did 
not subtract background from the data, because the timescale of the background variation 
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was much longer than 10^ s and did not affect the period determination of GX 1+4. The 
pulse period which gives the largest against a constant was searched. The time resolution 
of the data was 4 s in the 1990 observation, so we divided the folding pulse period into 28 
phase bins. The time resolution was 16 s in the 1991 observation and we divided the folding 
period into 7 phase bins. Finer binnings were not appropriate because the pulse period is 
about 10^ s. Heliocentric pulse periods we obtained are listed in Table 2. Average rates 
of pulse period change between the previous measurement and this measurement are also 
listed in Table 2. These results indicate that GX 1+4 was spinning down since 1987. 

EDITOR: PLAGE TABLE | HERE. 

The folded light curves are shown in Figure 1. The pulse profile of data set (a) is very- 
similar to those in 1988 and 1989: double peaked in low energy band and single peaked in 
high energy band accompanied by a sharp and deep dip. The sharp dip observed in data 
set (a) is not clear in data set (b). This may be due to the low time resolution of the data. 

EDITOR: PLAGE FIGURE H HERE. 

3.1.2. Spectral analysis 

GX 1+4 is located close to the galactic plane, hence it is difficult to accurately estimate 
the contribution from the galactic ridge emission (Koyama 1989 and the references therein). 
For this estimation, we sampled the spectra at two points in the sky, (/, b) = (3°.215, 

5°.744) and (5°.795, 7°.674); GX 1+4, (1°.943, 4°.785), is located on the line connecting 
these two points. The data were obtained as a part of galactic ridge observations, in which 
the data of GX 1+4 were also acquired. These two points are well separated from nearby 


bright sources such as the Kepler SNR and GX 9+9, and hence the spectra can be regarded 
as those of the pure galactic ridge emission. The spectra were htted by a model of thermal 
Bremsstrahlung and a Gaussian line with a center energy of 6.7 keV - the energy of Fe xxv 
Ka. The temperature of the thermal Bremsstrahlung was found to be 4.78 keV and 2.14 
keV, respectively, and the flux in unit of count rate at spectral peak was 2.5 ± 0.7 counts 
s“^ keV“^ and 0.95 counts s“^ keV“^, respectively. The Gaussian line was 0.89 ± 0.31 
counts s“^ at the former position and not detected at the latter. The flux of the galactic 
ridge emission at the position of GX 1+4 was then estimated by linear interpolation of 
the best-£t parameters. The resultant temperature and spectral peak flux were 6.1 keV 
and 3.27 counts s“^ keV“^, which corresponds to 1-37 keV flux of 18.1 counts s“^. The 
iron line flux was estimated to be 1.2 counts s“^. The ridge emission thus determined 
was subtracted from the energy spectra of GX 1+4 together with the intrinsic and Gosmic 
X-ray background utilizing the method described in Hayashida et al. (1989). The resultant, 
background-subtracted spectra are shown in Figure 2. 

EDITOR: PLAGE FIGURE g HERE. 

The energy spectrum of the observation in 1991 was well htted by a power law 
attenuated by a cold absorber, as was the case in the three previous observations in 1987, 
1988, and 1989. The spectrum of 1990 was, however, better expressed by a partial covering 
model of a power law continuum. A prominent iron line was also detected at about 6.4 keV. 
This line is considered to be a huorescence line of neutral or slightly ionized iron in the 
vicinity of the source, not a collisionally excited emission line of highly ionized iron whose 
center energy is at 6.7 or 6.9 keV. The huorescence line provides important information 
on the physics of the system (Makishima 1986). The center energy gives the ionization 
stage of the illuminated matter and thus the distance of the plasma from the source and 
the density. The geometry of the attenuating matter can be derived from the equivalent 


9 


width of the line and the attenuation column density. Therefore, it is essential to resolve 
a 6.4 keV line from 6.7 and 6.9 keV lines, and determine the central energy and equivalent 
width precisely. Unfortunately the 1991 energy spectrum has coarse energy binning and 
the line center energy could not be determined with enough accuracy to distinguish the 6.4 
keV and 6.7 keV lines. As seen later, ASCA data show that the iron line in GX 1+4 is 
centered at 6.4 keV. Thus, in the subsequent analysis of Ginga spectra, we hxed the line 
center energy to 6.4 keV. We summarize in Table 3 the best-£t parameters of both sets 
of observations. The spectra and the best-£t model are shown in hgure 2. We detected 
signihcant X-ray emission from GX 1+4 up to ~50 keV in data set 2. However, we could 
not hnd any signihcant structure of cyclotron resonant scattering in the energy spectrum. 

EDITOR: PLAGE TABLE | HERE. 

The absorption-corrected X-ray luminosity is calculated as L2_20 keV ~ ^ 

erg s“^ and 0.27 x 10^^ erg s“^ for 1990 and 1991, respectively, for the assumed distance of 
10 kpc. In 1989, the source was found to be as luminous as 10^^ erg s“h It was interpreted 
that GX 1+4 was recovering from the low state (Sakao et ah 1990). In 1991, however, it 
decreased again to the very faint level as in 1987, and was still in the steady spindown 
phase. 


3.2. ASCA data 

GX 1+4 was detected with both GIS and SIS with an average count rate of about 1 
counts s“^ (GIS: 1-10 keV). No signihcant time variation was seen except for the 122 s 
pulsation. Source photons were extracted from the circular region centered on the GIS 
image peak with a radius of 6 arcmin. As for SIS, photons from the whole chips (SOGl and 


S1C3) were used to get better statistics. Both timing and spectral analyses were done using 
these source photons. 


3.2.1. Timing analysis 

The pulse period of GX 1+4 was determined using the GIS data, which have superior 
time resolution. We applied neither detector-background subtraction nor dead time 
correction to the data, because such corrections are not signihcant for timing analysis of 
medium brightness sources such as GX 1+4. Using the epoch folding technique, we found 
the barycentric pulse period to be 122.019 ± 0.027 s. Folded pulse prohles are shown in 
hgure 3 in several different energy bands. Deep dips similar to that observed with Ginga are 
clearly seen in the pulse prohle. Below 3 keV, galactic ridge emission dominates and hence 
the pulse amplitude is much reduced. It should be noted that the dip structure in the pulse 
prohle becomes shallow in the energy band between 6.28-6.52 keV, where the huorescent 
iron line is dominant (see Figure 4). Iron line hux contributes about 40 % in this energy 
band and the reduction of the dip depth (~50 % in 6.28-6.52 keV compared to ~80 % in 
the other bands) is consistent with no pulsation in the iron emission line. This conforms 
to the scenario that the huorescent iron line originates from the spherically symmetric 
absorbing matter as discussed in section 4. The other energy bands besides <3 keV and 
6.28-6.52 keV show a similar pulse prohle and pulse fraction, indicating no pulse phase 
dependence of the spectral parameters, such as the photon index or the absorption column 
density. 


EDITOR: PLAGE FIGURE | HERE. 
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3.2.2. Spectral analysis 

We calculated a phase-averaged energy spectra of SIS and GIS using the source photons 
described previously. For the GIS spectrum, background photons were extracted using the 
GIS data from the region symmetric to the source region relative to the bore sight axis 
of the XRT. These background data include galactic ridge emission. We constructed an 
energy spectrum of the galactic ridge emission by subtracting the high latitude background 
data, which are publically available. We found that the background spectrum can be well 
represented by emission from an optically thin hot plasma model by Mewe et al. (1985). 
The best-£t temperature obtained was 6.1 keV with equivalent hydrogen column density of 
2.3 X 10^^ cm“^. Because we could not extract background data for SIS due to the small 
FOV, we just subtracted the high latitude blank sky data as background. Galactic ridge 
emission was included in a model spectrum for spectral fitting. 

We found that a single power law modified by cold matter absorption plus a narrow 
emission line can reproduce the energy spectrum. The model functions we employed are 

his{E) = fueUkT) exp{-aNi) + L E-^ + B exp 1 exp(-aiV^-) (1) 

\ -^width / _ 

and 

/gis(B) = \a E-^ + Bexp (-kLyETjl exp(-ffjvr) (2) 

for the SIS and the GIS data, respectively. Here, E is the X-ray energy, fueweikT) the model 
spectrum by Mewe et ah, F the photon index, a the photo-absorption cross section, iVn 
the hydrogen equivalent column density, Ec the line center energy, i^width the intrinsic line 
width, which is assumed to be effectively zero in the fitting, and A and B flux normalization 
factors. The model of the galactic ridge emission, fueweikT), is included only in the SIS, 
because it is subtracted from the GIS spectrum. The spectra and the best-fit model are 
shown in Figure 4, and the best-fit parameters are listed in Table 3. The luminosity of 
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GX 1+4 in 2-20 keV (same as the Ginga band) is estimated as 1.1 x 10^^ erg s ^ using the 
best-£t model for the assumed distance of 10 kpc. 

EDITOR: PLAGE FIGURE | HERE. 

To determine the ionization (state) of iron, we htted the 5-10 keV SIS data with a 
power law attenuated by the iron K edge plus iron Ka and K/3 lines. The line-flux ratio of 
the iron K/1 to Ka was hxed to the value expected in fluorescence of Fe l. The center-energy 
ratio of the Kj3 to Ka was also fixed to the value expected in Fe l. We tested several values 
of this ratio, and found that they do not change the result. For SISO and SISl, the center 
energy of the iron Ka line was htted to be 6.400]'^o;oJo keV and 6.440lo;o}2 keV, respectively, 
and the energy of the iron K edge was htted to be 7.163lo;o72 keV and 7.216lo;o72 keV, 
respectively. The errors are 90 % conhdence intervals. The line center energies of both 
sensors were inconsistent with each other. This is considered to be due to gain-calibration 
uncertainties of the SISs caused by the decreasing charge transfer efficiency (e.g., Gendreau 
1995; Yamashita et ah 1997). For measurement of the ionization degree free from the gain 
uncertainties, we utilize the ratio of the edge energy to the line center energy. The ratios 
were 1.119lo;oii and 1.120lo;oii for SISO and SISl, respectively. It should be noted that 
these values are consistent with each other within the statistical errors. Thus the averaged 
energy ratio was accurately determined to be 1.120i'io;oo8, which corresponds to Fe l-Fe iv 
(Lotz 1968). 


3.2.3. Image analysis 

We searched for a dust scattering halo around GX 1+4. X-rays from a source near the 
galactic plane such as GX 1+4 would suher from scattering by interstellar matter and may 
produce a halo around the source. Efficiency of the dust scattering is a strong function of 
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X-ray energy and the dust scattering halo is most prominent at lower energies. To search 
for the halo, we compared the radial prohle of the image above and below 5 keV. The 
dividing energy was selected to be as low as possible while retaining statistical signihcance 
in the lower energy band. This method is unaffected by possible attitude fluctuation of the 
satellite during the observations. Although GX 1-1-4 is located very close to the direction 
of the galactic center, the effect of dust scattering was not seen in the image. This is not 
surprising, because the source is heavily absorbed and has little X-ray emission in the lower 
energy band, where the dust scattering is most prominent. Predehl, Friedrich, & Staubert 
(1995) found an X-ray halo around GX l-|-4in a ROSAT-PSPG image below 2.4 keV. The 
halo extended to a few arcmin and was brighter than GX 1-1-4. Predehl et ah suggested 
that this halo was the reflection of the flare 70 hours before the ROSAT observation. That 
is consistent with our result. 

4. Discussion 

Through the Ginga and ASGA observations, we found large changes in the energy 
spectrum, especially for the absorption column density and the iron emission line. On 
the other hand, the spindown rate decreased slightly in 1990 and 1991 followed by an 
increased rate, and GX 1-1-4 entered a short episode of spinup in 1994 just after the ASGA 
observation. The spinup episode is recognized as a glitch in 1995 in Figure 1 in Ghakrabarty 
et al. (1997), or a flare beginning at MJD 49638 in Figure 2 ibid. In this section, we consider 
the possible connection between the circumstellar matter and the spin period change. 
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4.1. Geometry of the circumstellar mater 

Hydrogen column densities obtained from the GX 1+4 observations with Ginga and 
ASGA range between 10^^ — 10^^ cm“^. These column densities are several orders of 
magnitude larger than the galactic absorption of 6 x 10^^ cm“^ (Dickey & Lockman 1990), 
and are considered to have a circumstellar origin. The equivalent width of the fluorescent 
iron line is useful to deduce the distribution of the circumstellar matter. We show in 
Figure 5 a scatter plot of the equivalent width of the iron line and the hydrogen column 
density obtained from the ASGA and all the Ginga observations. The relation expected 
from a spherical distribution of matter is also indicated in the figure. The column density 
drastically increased to 10^“^ cm“^ in 1991 with the increase of iron equivalent width to 3 
keV. This equivalent width is very close to that expected from a spherical distribution of 
the matter. The ASGA observation showed that the column decreased in 1994 to a level 
comparable to that in 1990. When the column density is relatively low, the data points are 
scattered around the relation expected from the spherical distribution; but, the deviation is 
at most factors of a few. This scatter may be due to a clumpy structure of the circumstellar 
matter. In fact, we found that a partial covering model gives a better fit to the energy 
spectrum in 1990. This supports the assumption that the circumstellar matter has a clumpy 
distribution. Except for this dumpiness, we can consider that the circumstellar matter has 
a spherical symmetric distribution around GX 1+4. 

EDITOR: PLAGE FIGURE | HERE. 

We found from the ASGA observation that the ionization degree of iron in the 
circumstellar matter is very low (Fe l-Fe iv). Thus, the ^-parameter (= Lx/nr^), which 
determines the photo-ionization degree, may be smaller than 30 erg cm s“^. With 
the knowledge of Lx and Nh (~ nr), we can estimate the characteristic scale of the 
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attenuating-matter distribution or the typical distance of the matter from the source. 
ASCA data are not very useful to evaluate Lx, because the data cover only below 10 
keV. According to balloon experiments, the hard X-ray luminosity (20-100 keV) was 
approximately ~ 2 x 10^^ erg s“^ (for the distance of 10 kpc) during the 1990’s except 
for the bright state in 1993 December (Rao et ah 1994; Paul et ah 1997). Therefore, we 
assume Lx = 2 x 10^^ erg s“^ as the total X-ray luminosity during the ASCA observation. 
This X-ray luminosity and the column density of iVn = 2.2 x 10^^ cm“^ give the radius of 
the attenuating circumstellar matter to be larger than 3 x 10^^ cm. Thus, it is concluded 
that the source was enveloped in a spherical gas cloud with a radius of > 3 X 10^^ cm. 
Considering various uncertainties in the above calculation, the result should be regarded 
as an order of magnitude estimation. However, it is still larger than the co-rotation radius 
Tco = 4 X 10® cm and the Alfven radius ta = 4 x 10® cm even for the extraordinarily strong 
magnetic field of 10^"^ G. Instead, it may be comparable to the size of the Roche lobe, if the 
binary period is of order a year. 


4.2. X-ray heating of the circumstellar matter 

GX 1-1-4 is thought to accrete matter from the slow stellar wind of the companion. 
The large size of the circumstellar matter and its nearly spherical distribution around 
GX 1-1-4 is naturally interpreted as the dense stellar wind from the companion captured in 
the Roche lobe of GX 1-1-4. The temperature of the plasma is estimated as a few eV from 
the value of the .^-parameter (~ 30 erg s“^ cm), which corresponds to a thermal velocity 
of about 30 km s“^. On the other hand, the Kepler velocity at 10^^ cm from the neutron 
star is about 100 km s“^. Thus the circumstellar matter is considered to be bound to 
GX 1-1-4. It is noteworthy that, if GX 1-1-4 becomes slightly brighter, the circumstellar 
matter can easily escape from GX 1-1-4. Suppose GX 1-1-4 becomes 3 times more luminous. 
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the .^-parameter of the circumstellar matter would increase up to 10^. This corresponds 
to a plasma temperature of ~30 eV, and hence a thermal velocity of 100 km s“^. This is 
comparable to the Keplerian velocity and the plasma can escape from GX 1+4. 

During the spin down phase of GX 1+4 since 1985, the stellar wind captured by 
GX 1+4 may be stagnated around the neutron star at a distance of 10^^ cm or more. The 
low luminosity characteristic of the spin down phase would reduce the amount of matter 
escaping from GX 1+4 by X-ray heating, resulting in a larger column density. Some other 
effects such as a change of the wind velocity may also have contributed to the increase 
in column density. Because the amount of the stellar wind captured by the neutron star 
depends on the wind velocity as oc a small change of the wind velocity can result in 
the large change of the captured mass. 

When GX 1+4 spins up, as in 1970’s and in the short episode of 1994-1995, the X-ray 
luminosity is at least several times larger than that in the spindown phase. The X-ray 
luminosity may be large enough to heat up the circumstellar matter and to expel much of 
it from the vicinity of GX 1+4. This explains the relatively low column density observed 
during the 1970’s (Becker et ah 1976). 

Just before the short episode of spinup in 1994-1995, the column density was observed 
with ASGA to be reduced to the level of 1990. This means that most of the accumulated 
circumstellar matter had been expelled from the vicinity of GX 1+4 before the transition 
to spinup. Thus it is suspected that the accumulation of the circumstellar matter may not 
be the cause of the spinup, and is just controlled by the X-ray luminosity. This may be 
understood if the matter which affects the spin period change is located near the Alfven 
radius, not at the distance of 10^^ cm from the neutron star. 
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4.3. Optical observations 

Optical observation is powerful for plasma diagnosis, and it is important to check 
whether our interpretation of the X-ray data is consistent with optical data or not. Here we 
briefly describe a preliminary analysis of optical observations, and examine the consistency 
between the two bands. We do not discuss in detail as this will be reported in a forthcoming 
paper by Greenhill (1998). 

V2116 Oph, the optical counterpart of GX 1-1-4, was observed with the 3.9 meter Anglo 
Australian Telescope and the RGO spectrograph on 1994 September 25, 11 days after the 

O 

ASGA observations. The measurements covered the wavelength range 3700 to 8900 A with 

O 

a resolution of 3.2 A. On September 26 measurements centered on the Ha wavelength were 

O 

made with a resolution of 1.0 A FWHM. Photometric observations in the V, Ha and I 
bands were made with the Mount John University Observatory 1 meter telescope on hve 
nights (August 25 and 27, September 2, and October 3 and 5) covering the time of the 
ASGA observations. The Ha flux increased slowly during this period and was about 0.4 
magnitudes brighter on October 5. The V-I color changed only marginally. In view of these 
relatively small changes we believe our AAT data are representative of conditions in the 
source during the ASGA observations. 

The measured and de-reddened fluxes of some prominent emission lines are listed in 
Table 4. The reddening correction assumes the same interstellar extinction as observed 
by Ghakrabarty and Roche (1997). The observed Ha flux (1.6 x 10“^^ erg cm“^ s“^) was 
about half the value observed by Davidsen et ah (1977) during the high state in the mid 
1970’s and by Ghakrabarty & Roche (1997) in 1993. The measured Ha/H/5 line strength 
ratio of ~30 was less than 25 % of the ratio measured by Davidsen et ah during the high 
state. Numerous Balmer, Paschen, He l, O l, O ll, [O ill], Na l,Ga ll and Fe ll lines were 
seen in emission but higher excitation lines were weak or undetected. He ll 4686 was weak. 
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the C II /N III feature common in X-ray sources at 4640-4650 A was not detected; nor was 

O 

the 6830 A Raman line although this was very strong during the high state and in 1988 
(Davidsen et ah 1977; Chakrabarty & Roche, 1997). The O l 8446 line, an indicator of Ha 
optical depth (Netzer & Penston 1976) was almost an order of magnitude fainter relative to 
Ha than during the 1993 measurements by Chakrabarty & Roche (1997). 

EDITOR: PLACE TABLE | HERE. 

Following Chakrabarty & Roche (1997) we assume that the emission lines are excited by 
photo-electric interactions in circumstellar matter by UV emission from the disk. Further we 
adopt the hypothesis developed earlier that the matter is trapped in the gravitational held 

O 

of the neutron star. The Ha line width (~2 A FWHM after correction for the spectrograph 
resolution) is consistent with Doppler broadening due to matter gravitationally bound to 
the neutron star at a radius of ~ 3 x 10^^ cm. As noted by Chakrabarty & Roche (1997) 
the presence of Fe ll and absence of Fe ill lines implies a temperature in the emission region 
of <30,000 K. The ratio of the corrected O l 8446 to Ha emission is ~0.01 indicating an 
Ha optical depth ~150 (Netzer & Penston 1976). From this and the corrected line strength 
ratios of the hrst 4 Balmer lines we hnd, using the model of Drake & Ulrich (1980), an 
electron density Ug ~ 3 x 10^° to 10^^ cm“^ and a temperature ~20,000 K in the emission 
line region. Similar values for the electron density and plasma temperature are obtained 

O 

from the ratios of the He l (5876, 6678 and 7065 A) emission lines (Proga et ah 1994). 
These values are consistent with the hydrogen equivalent column density A^h = 2 x 10^^ 
cm“^ and plasma temperature derived from the ASCA spectroscopy assuming a spherical 
distribution of matter with a radius of ~ 3 x 10^^ cm about the source. The densities are 
one to two orders of magnitude greater than those estimated by Davidsen et ah (1977) and 
by Chakrabarty & Roche (1997) and the radius is an order of magnitude smaller. Thus it is 
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shown that the parameters of the system derived from the X-ray data, assuming a spherical 
distribution of circumstellar matter coincide with those of the optically emitting plasma. 

Since the separation of the neutron star and its giant companion is expected to be 
greater than 1 AU (1.5 x 10^^ cm) the matter in our model can be bound to the neutron 
star. Clearly, in the Davidsen et al. (1977) model the emission region surrounds the whole 
binary system. 


4.4. A negative feedback mechanism 

The Ginga and ASCA observations suggest that gas is accumulated in the vicinity of 
the source, probably by the stellar wind from the companion star, and may be blown off by 
an increase in X-ray luminosity. If the accreting matter is supplied from the stellar wind as 
suggested (e.g., Dotani et al. 1989; Davidsen et al. 1977), even a small increase in X-ray 
luminosity would lead to a decrease in accretion. This constitutes a negative feedback 
process and suppresses the X-ray luminosity. The situation is similar to the self-regulation 
of the luminosity at the Eddington limit. GX 1-1-4 entered a low state in the early 1980’s, 
and, except for occasional short flares (Chakrabarty et al. 1997), remained faint until the 
ASCA observation in 1994. During the spindown phase, the activity may be suppressed 
to a low level by the negative feedback mechanism. The condition for the mechanism to 
work is rather difficult to satisfy compared to Eddington accretion; (1) The .^-parameter 
of the accretion flow needs to be high enough so that the thermal velocity can exceed 
the escape velocity. In the case of GX 1-1-4, the relatively low density of the accumulated 
matter, > 10^^ cm“^ in 1991, resulted in a large .^-parameter. On the other hand, in binaries 
where the accretion matter is supplied by Roche-lobe overflow, the density ~ 10^® cm“^ is 
presumably too high for the mechanism to work. (2) The flow velocity must be lower than 
the thermal velocity attainable through X-ray heating. If the ram pressure of the flow is 
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higher than the thermal pressure, trapped matter will not be blown off by X-ray heating. 
The typical wind velocity of M giants is 10-30 km s“^ and it is comparable or smaller than 
the thermal velocity of a plasma of ~ 30. When GX 1-1-4 is in a flare state and ^ > 100, 
the thermal velocity exceeds the wind velocity. In the case of massive binaries with early 
type companion stars such as Cen X-3 or Vela X-1, the wind is supersonic and the negative 
feedback mechanism is not likely to work, even if the sources are wind fed. 

Chakrabarty, van Kerkwijk, & Larkin (1998) reported that the wind velocity of 
GX 1-1-4 is as large as 250 ± 50 km s“^. Ghakrabarty et al. suggested that radiation from 
the accretion disk or the neutron star may contribute to the acceleration of the outflow. 
That is consistent with the picture of the negative feedback, which predicts an outflow of 
circumstellar matter heated by the source. However, if the initial velocity of the stellar 
wind were as large as 250 km s“^, our model would fail. 

It should be noted that this negative feedback will not keep the X-ray luminosity 
constant. The X-ray source blows off only matter located far from the source with small 
gravitational binding energy. Therefore, the fuel supply decreases after a certain time 
comparable to the orbital period. Assuming the central mass to be 1.4 Mq, the orbital 
period is estimated to be 10^ days at the distance of 10^^ cm. Such a negative feedback 
with a time delay would result in temporal variability such as oscillations or limit cycles. 
This mechanism predicts that; large X-ray flares cannot last long because they even blow 
off matter close to the source, and that; small flares or weak activities last longer. That is 
qualitatively consistent with the pulsed flux history reported by Ghakrabarty et al. (1997). 
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4.5. Spindown and spinup phase 

Our interpretation of the variability in the luminosity and the absorbing column density 
through the negative feedback mechanism explains the very steady spindown trend. The 
mass supply, averaged over a time longer than the feedback timescale, would not change 
very much even if the mass loss rate of the companion changes considerably. Either the 
retrograde disk or a propeller effect processes can spindown the neutron star steadily under 
the feedback mechanism. 

In the spinup phase during the 1970’s and in 1994, the feedback mechanism obviously 
did not work. In the 1970’s, the luminosity was typically as high as 10^^ — 10^® erg s“^ 
and was enough to blow off a stellar wind with a velocity of 10 km s“^ and a density of 
< 10^^ cm“^. The wind could not reach the surface of the neutron star unless it was as 
dense as > 10^^ cm“^. To avoid blow off, the accreting matter had to be supplied by a 
very dense stellar wind or by a Roche-lobe overflow. The condition of transition between 
spindown and spinup phases is not understood yet. Since the companion is a low mass star. 
X-ray irradiation of the surface may affect the thermal structure and the evolution (e.g., 
Podsiadlowski 1991; D’Antona & Ergma 1993). Harpaz & Rappaport (1994) show that 
such a star may repeat short episodes (1-100 yr) of mass transfer. It is possible that the 
spinup and spindown phases of GX 1-1-4 are part of the episodic mass transfers induced by 
X-ray irradiation. 

5. Conclusion 

GX 1-1-4 in spindown phase was observed with Ginga in 1990 and 1991, and with 
ASGA in 1994, a month before the transition to a short spinup episode. The behavior of 
the source in 1990 and 1991 can be summarized as below: 



22 


• It was found that the absorption column density drastically increased to 
A^h = (1-4 ± 0.2) X 10^^ cm“^ in 1991, indicating a rapid accumulation of 
matter in the vicinity of the source. 

• The luminosity decreased to T 2 _ 2 q = 0.72 x 10^'^ erg s“^ in 1990, and to 

0.27 X 10^^ erg s“^ in 1991, which is almost the same level as the lowest activity in 
1987. Nevertheless, the spindown rate was very stable at (3 — 4) x 10“® s s“^ during 
the time. 

Properties of the source deduced from the ASCA observation can be summarized as below: 

• The absorption column density decreased to Ah = (2.08 ± 0.02) x 10^^ cm“^ in 1994 
after a long interval in a high absorption phase. 

• The ionization stage of iron in the absorbing matter was derived to be Fe l-Fe IV from 
the ratio of the absorption-edge energy to the emission-line energy. 

• The .^-parameter of the absorbing matter was <30, consistent with a spherical 
distribution with a radius of > 3 x 10^^ cm about the source. 

• The luminosity brightened to T2_20 keV “ ^ which is the brightest in 

the spindown phase, and almost comparable to that in the highly-active spinup phase. 

The picture of GX 1-1-4 in a gas cloud is also supported by the optical observations covering 
the ASCA observation. The Doppler broadening of the Ha line was consistent with the 
Kepler motion at ~ 3 x 10^^ cm from a neutron star, and electron density and temperature 
estimated from the line strength ratios were consistent with the values expected from the 
picture. 

To explain the accumulation and the disappearance of the absorbing matter, a negative 
feedback mechanism stabilizing the spindown rate is suggested. Gas supplied by a slow 
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stellar wind from the companion star would be blown off, if the luminosity of the neutron 
star is high. When the luminosity is low, the neutron star would accrete and increase the 
luminosity. Therefore, the accretion rate and thus the spindown rate may be stabilized 
under a certain level. This mechanism does not depend on the nature of the unknown 
spindown mechanism, for which a retrograde disk and high magnetic held were proposed. 

Since this negative feedback mechanism will not explain transitions between spinup and 
spindown, the transitions should be explained by some other mechanism such as episodic 
mass transfers. 

TK would like to thank Dr. T. Aoki, Dr. K. Asai and Dr. H. Negoro, for helping to 
analyze Ginga data. TK is supported by the Special Postdoctoral Researchers Program of 
RIKEN. 
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Fig. 1.— Folded light curves of the Ginga data in (a) 1990 September (1-37 keV) and (b) 
1991 September (2-60 keV). Background and galactic-ridge components were not subtracted. 
See Table 2 for the heliocentric pulse period of these data. 

Fig. 2.— Energy spectra taken with Ginga in (a) 1990 September and (b) 1991 September. 
The intrinsic and the galactic backgronnd are snbtracted. The best-ht model (solid) and the 
contribution of the power law (dot) and the iron line (dot-dash) are plotted as histograms. 
The contribntion of leaked power law (dash) is also plotted in (a). See Table 3 for the best-ht 
parameters. 

Fig. 3.— Folded light curves of the SIS data in six energy bands in 1994 September. The 
temporal average of the hnx of each curve is normalized to nnity. Both sensors are combined. 
The epoch at phase 0 is 49610.00 (MJD). Backgronnd and galactic-ridge components were 
not snbtracted. 

Fig. 4.— Spectra taken with ASGA (a) SIS and (b) GIS. The best ht model is plotted 
as histogram. See Table 3 for the best ht parameters. From the GIS data, the galactic 
ridge emission is snbtracted, while only the dihuse cosmic X-ray backgronnd component is 
subtracted from the SIS data. 

Fig. 5.— History of the equivalent width and the absorption colnmn density. The data before 
1989 were taken from Sakao et al. (1990). The colnmn density of 1990 was mnltiplied by the 
covering fraction. The relation of eqnivalent width and absorption colnmn density expected 
in the case of isotropic distribntion of absorbing matter (Makishima 1986) is plotted as a 
dashed line. 
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Table 1: Journal of observations 


No 

Satellite 

Start Time 

(UT) 

End Time 

(UT) 

Exposure 

(ks) 

Energy Range 

(keV) 

1 

Ginga 

1990/9/4 13:47 

1990/9/5 04:28 

10.8 

1-37 

2 

Ginga 

1991/9/14 16:56 

1991/9/16 01:41 

37.8 

2-60 

3 

ASGA 

1994/9/14 23:53 

1994/9/15 22:00 

30.0 

0.5-10 







Table 2: Results of timing analysis 


No 

Mission 

Epoch 

(MJD) 

Pulse period 

(s) 

Period derivativef 

(10-® s s-i) 

1 

Ginga 

48138.0 

114.63 ±0.02| 

3.1 ±0.1 

2 

Ginga 

48513.0 

116.17 ± 0.031 

4.7 ±0.2 

3 

ASGA 

49610.5 

122.02 ±0.03§ 

6.17 ±0.03 


fAverage between observations. 
IHeliocentric value. 
§Barycentric value. 








Table 3: Spectral parameters of GX 1+4 


No 

Power Law 

Absorption 

Iron Line 


Normalization^ 

Photon index 

Nh 

EW 

Genter 




(10^^ cm“^) 

(keV) 

(keV) 

1 

0.156 + 0.026 

1.98 + 0.07 

3.0+ 0.3 

0.11 + 0.03 

6.4 (fixed) 




At=0.87 



2 

0.022 + 0.015 

1.49 + 0.23 

14 + 2 

3.1 + 0.2 

6.4 (hxed) 

3 

0.0503 + 0.00261 

1.19 + 0.09t 

2.15 + 0.03t 

0.094+ 0.005t 

6.42 + 0.04§ 


lln unit of photons cm ^ s ^ keV ^ at 1 keV 
jPartial covering factor. 

I Average between SIS and GIS. 

§SIS only. Including systematic errors. 








Table 4: Prominent emission lines 


Ion 

Wavelength 

O 

(A) 

Measured Flux 

(10“^® erg cm“^ s“^) 

Corrected Fluxf 

(10“^"^ erg cm“^ s“^) 

H I 

6563 

1590 

619 

H I 

4861 

52 

122 

H I 

4341 

6.6 

31.4 

H I 

4102 

2.6 

17.2 

He I 

5876 

39 

28.3 

He I 

6675 

27 

9.6 

He I 

7065 

54 

14.4 

He II 

4686 

3.2 

9.5 

0 I 

8446 

50 

6.0 


Measured at the AAT on 1994 September 26. 


fThe reddening corrections are those used by Chakrabarty and Roche 1997 
















































































